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Abstract 

The BKL conjecture, stated in the 60s and early 70s by Belinski, Khalatnikov and Lifshitz, 
proposes a detailed description of the generic asymptotic dynamics of spacetimes as they approach 
a spacelike singularity. It predicts complicated chaotic behaviour in the generic case, but simpler 
non-chaotic one in cases with symmetry assumptions or certain kinds of matter fields. 

Here we construct a new class of four-dimensional vacuum spacetimes containing spacelike 
singularities which show non-chaotic behaviour. In contrast with previous constructions, no 
symmetry assumptions are made. Rather, the metric is decomposed in Iwasawa variables and 
conditions on the asymptotic evolution of some of them are imposed. The constructed solutions 
contain five free functions of all space coordinates, two of which are constrained by inequalities. 

We investigate continuous and discrete isometries and compare the solutions to previous 
constructions. Finally, we give the asymptotic behaviour of the metric components and curvature. 
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1. Introduction 

1.1. Singularities in general relativity 

When Albert Einstein presented his theory of General Relativity in 1915 he did not give any 
non-trivial exact solutions to its field equations. Due to the complicated non-linear structure of 
the equations, he did not expect any to exist and calculated physical predictions using pertur¬ 
bation theory [T]. To his surprise, less than a month later, Karl Schwarzschild sent him a letter 
containing the Schwarzschild metric, a spherically symmetric solution of the vacuum Einstein 
equations. It is given, in Schwarzschild coordinates, as 

1 _ 2^^ ' dr^ + r^idO^ + sin^ 9 d^) • 

This solution contains, in these coordinates, an apparent singularity at r = 2m where the rr 
component of the metric diverges. This is the event horizon, which Schwarzschild set as the origin 
of his coordinate system. The solution also contains a real singularity at r = 0 which was found by 
David Hilbert in 1917. Hilbert considered both singularities real, as they could not be removed by 
an everywhere smooth and invertible coordinate transformation. In hindsight his requirement was 
too strict: The fact that in Schwarzschild coordinates the rr component of the metric diverges 
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at r = 2m simply means that these coordinates are badly chosen, indeed to transform from 
coordinates which don’t show the apparent singularity to Schwarzschild coordinates requires a 
transformation which diverges at r = 2m. 

In 1921 and 1922 Paul Painleve and Allvar Gullstrand independently discovered a spherically 
symmetric vacuum solution containing only a single singularity at r = 0 [113]. It was not realized 
at the time that this solution can be obtained from the Schwarzschild one by a coordinate 
transformation, i.e. it describes the same physical spacetime. This was finally discovered by 
Georges Lemaitre in 1932, who also correctly identified the r = 2m singularity as an apparent 
singularity caused by the choice of coordinates [Ij. The singularity at r = 0 cannot be removed 
by a coordinate transformation as the Kretschmann scalar, given by diverges there. 

This is a scalar quantity, constructed by contracting all indices of the Riemann tensor with itself, 
and is therefore independent of the chosen coordinate system. 

Despite this advance, the status of real singularities, such as the one appearing in the 
Schwarzschild or the cosmological FLRW solutions, was unclear. It was widely believed that 
they were an artifact of the symmetry assumptions made to obtain explicit solutions and had no 
relevance for the real world |5|. The idea was that, similarly to the Newtonian case, if matter was 
not perfectly symmetrically rushing towards a central point, the resulting angular momentum 
would prevent the formation of a singularity. 

The singularity theorems of Penrose and Hawking mm proved the opposite. They state 
that, given a trapped surface, an energy condition, and an assumption on the global structure 
of spacetime (e.g. no closed timelike curves), a singularity, in the sense of geodesic incomplete¬ 
ness, has to form. As small perturbations of an explicit solution containing a singularity would 
preserve the trapped surface, the perturbed solution also contains a singularity. These theorems, 
however, do not give any information about the nature of the predicted singularities, or about 
the behaviour of the metric near them. Indeed they do not even predict diverging curvature, 
only the existence of some geodesics, which cannot be extended beyond a finite value of the affine 
parameter along them. 

1.2. The BKL conjecture 

In a series of works, beginning in 1963, Belinski, Khalatnikov and Lifschitz (BKL) conjectured, 
based on heuristic arguments, that the dynamics of a generic spacetime containing a spacelike 
singularity would drastically simplify when the singularity is approached mm- They claimed that 
time derivatives of the metric would dominate compared to space derivatives, causing different 
spatial points to effectively decouple and turning the Einstein equations into a system of ODEs 
at each point. The solution of these ODEs is a generalisation of the Kasner metric, an explicit, 
homogeneous (but anisotropic) solution of the Einstein equations describing a spacetime which 
expands in some directions and contracts in others. It is given by 


d 



( 1 . 1 ) 


with the constants pj fulfilling pj = 1 and J2j p] = ^ (these conditions imply that at least 
one of the pj has to be negative, unless one is 1 and all others 0). The Kasner metric contains a 
singularity at t = 0 as the Kretschmann scalar, given by 


K = R'^^^^R^p^s = 16 {P3^-P3^) , 


diverges there. 
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The behaviour predicted by BKL consists of a series of time periods (often referred to as 
Kasner epochs) during which the metric behaves at each spatial point as the Kasner metric, 
but with spatially varying exponents. At the end of a Kasner epoch the Kasner exponents pj 
change rapidly to a new configuration causing an “oscillation” as previously expanding directions 
contract. As the singularity is approached, the Kasner epochs get shorter and shorter and the 
transition between epochs becomes sharper. 

Chitre [TUI Misner [TT] introduced a representation of the BKL behaviour as a (chaotic) 
billiard motion in an auxiliary space of the same number of dimensions as the space part of the 
spacetime. A “particle”, representing some parts of the metric, moves along straight, null, lines in 
a flat Lorentzian space and is elastically reflected off of (asymptotically) inhnitely high potential 
walls. The straight line motion represents a Kasner epoch while the (asymptotically) sharp 
reflections correspond to the transitions between epochs. This billiard approach is described in 
detail by Damour, Henneaux and Nicolai in [12] . 

Rigorous results concerning this chaotic case of the BKL conjecture are sparse: The only 
known example of a spacetime which shows the full chaotic BKL behaviour was constructed by 
Berger and Moncrief [TU]. They applied a solution generating transformation to a homogeneous 
cosmological solution, yielding a 17(1) symmetric one. The resulting solution shows chaotic be¬ 
haviour but it is very restricted, containing no free functions, and only three arbitrary constants. 

Numerical investigations do, however, provide strong evidence supporting the BKL conjecture 
M- More recent simulations have shown that, while generically the spatial derivatives do become 
negligible, there are exceptional points at which they instead increase exponentially, giving spikes 
in the metric components m- In the class of Gowdy spacetimes explicit (non-chaotic) solutions 
exhibiting this behaviour have been found m- The appearance of these spikes hints at more 
complicated detailed behaviour within the general dynamics predicted by BKL. 

1.3. Asymptotically simple behaviour 

Belinski and Khalatnikov argued that coupling a massless scalar field to the Einstein equations 
would reduce the BKL behaviour to a simpler, non-oscillatory one, described by a single Kasner 
epoch, which is sometimes called AVTD (Asymptotically Velocity Term Dominated) [T7]. This 
was rigorously proven, including the case of a stiff fluid, by Andersson and Rendall [TU] . 

If a p-form held is added to the scalar one, the resulting behaviour is either simple (single 
Kasner epoch) or chaotic, depending on the coupling constant between them. This was shown 
by Damour, Henneaux, Rendall and Weaver [TU] . 

In the billiard picture, the addition of matter increases the dimension of the auxiliary space, 
as the particle describes not only the metric components but also the values of the matter helds. 
In addition, the evolution equations for the matter helds add additional potential walls. If a 
null line in the auxiliary space, which does not intersect any of the walls, exists, the resulting 
behaviour is simple, as a single Kasner epoch lasts up to the singularity. 

The addition of matter is not necessary for non-chaotic behaviour: Demaret, Henneaux and 
Spindel [20j argued, using similar heuristic arguments as BKL, that in 10 or more spatial dimen¬ 
sion AVTD behaviour is generic. 

Even in lower dimensions, where BKL predict chaotic behaviour in the generic case, solutions 
which show non-chaotic behaviour exist. They are characterized by symmetry assumptions or 
conditions on their asymptotics. These assumptions cause some of the potential walls in the 
billiard picture to vanish at least asymptotically. 

This reduction was first proven for the polarized Gowdy subclass of the symmetric space- 
times by Chrusciel, Isenberg and Moncrief EH Eg. It was later extended to a larger class 
of Gowdy spacetimes by Kichenassamy and Rendall [23] using a newly introduced “Fuchsian” 
method, which was then applied to more general symmetric spacetimes by Isenberg and 
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Kichenassamy [ 21 ]. An extension to so-called “half-polarized” spacetimes was achieved by 
Clausen and Isenberg [23 ■ Ames, Beyer, Isenberg and LeFloch [26] extended the previous results 
on spacetimes to lower regularity. All results on symmetric spacetimes focused on the 
case of 3 -I- 1 dimensions. 

17(1) symmetric AVTD solutions, with only one Killing field, were constructed by Isenberg, 
Moncrief and Choquet-Bruhat in 3 -I- 1 dimensions [ 2711251111 ]. 

The results obtained using Fuchsian methods do not necessarily provide generic solutions in 
the class of metrics considered, only the existence of families of solutions which contain a number 
of arbitrary functions. As these functions specify the asymptotic behaviour of the solution, there 
is no obvious link with functions in the initial data. Within the class of Gowdy spacetimes, 
genericity of AVTD behaviour was proven by Ringstrom |30j . 


1 . 4 . This work 

All previous results on simple behaviour in the vacuum case were obtained by starting with 
an ansatz for the metric which included one or more continuous symmetries. In the billiard 
picture this causes one or more of the walls to vanish identically at all times. 

Here a new class of non-chaotic vacuum solutions will be constructed without starting from 
such an ansatz. Instead, the decay of certain parts of the metric, defined by writing it in so- 
called Iwasawa variables, will be required. This causes some of the walls in the billiard picture 
to vanish asymptotically. The approach is based on work by Damour and de Buyl m who gave 
a precise statement of the BKL conjecture using this decomposition of the metric. Their work is 
an extension of m by Damour, Henneaux and Nicolai. The new class of solutions includes the 
polarized Gowdy ones, but not the other classes mentioned above. It is at the same time more 
general, as it includes free functions which depend on all space coordinates, and more specific, 
as some asymptotically free functions in e.g. the “half-polarized” case are here assumed to 
become constants in space. 

In sections!^ to the relevant parts of [3T] are described: Section [^describes the conventions 
and choice of gauge used and introduces the Iwasawa decomposition of the metric. In sections]^ 
and [^ the action and Hamiltonian in the Iwasawa variable form are given, including the potential 
“walls”. Section [^ states the Fuchs theorem, which is the main tool used in the construction 
of the new class of solutions. In section [7] the evolution equations are written in Iwasawa form 
and the approach to constructing solutions with specified asymptotic behaviour, as used e.g. by 
Rendall, is detailed. Finally, in section]^ the new class of solutions is constructed. In section]^ 
the new solutions are analysed: In sections |9.I| and |9.2| possible isometries of the solutions are 
investigated and in section [9.3| their relationship with previously known classes is described. 


In the appendices some of the calculations are given in more detail: In Appendix A the deriva¬ 


tion of the Iwasawa form of the Hamiltonian is given in full. [Appendix B 


contains comments on 


the form of the evolution equations used. In [Appendix G the Iwasawa form of the momentum 
constraint equations is derived, following m- In [Appendix D[ the evolution equations for the 


constraints are derived in the chosen gauge. [Appendix E gives the asymptotic behaviour of the 
metric components and curvature for the new class of solutions. [Appendix F 


shows that the 


new solutions can be constructed for arbitrary values of the cosmological constant. 
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2. Conventions, Iwasawa decomposition 

We work with a signature. Greek indices .. run from 0 to D = d + 1, 

Latin ones a,b,c,... from 0 to d. The metric (in D = d + 1 dimensions) is written in the form 

ds^ = —N{t, a;*)^dT^ + (r, x^)dx'‘dx^ , 

i.e. with vanishing shift vector and lapse N{t,x'^) = y^det gij. This pseudo-gaussian gauge 
has the unusual property that changes of the spatial coordinates also change the slicing of the 
spacetime. 

The spatial metric is then decomposed into Iwasawa variables j3°' and Af® as 

a—1 

Here the /3“ and are functions of all coordinates (including time) and A/”® vanishes for 
all a > z and is 1 for a = z (i.e. Af^ is upper triangular with ones on the diagonal). As the 
determinant of Af is 1, the determinant of the spatial metric only depends on the /3“ and is given 
by 

det5 = e“^^“^ . (2.1) 

The /3“ are referred to as “diagonal degrees of freedom” while the A/”® are the “off-diagonal 
degrees of freedom” (in fact both are relevant for all the metric components except gii). 

This decomposition corresponds to a Gram-Schmidt orthogonalization of the coordinate 
coframe dx*. 

The Iwasawa variables j3^ and Af^ have the advantage that they explicitly separate parts 
of the metric which have different asymptotic behaviour: As we will see later, the Af® go to 
constants as r —)■ oo while the /3“ approach linear functions. 

The Iwasawa coframe and its dual are defined as 

6(“=Af“dx* and ea = {N-^)\d^. (2.2) 

The structure functions of the Iwasawa coframe, denoted , are defined by 

dr = - ^ cV ^ [66, ee] = ea , 

and are given in terms of the A/”® as 

GV = E • (2-3) 

i,k 

In the 6 ^ coframe the metric takes the diagonal form 

d 

gydx^dx^ = E r . 

a—1 
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3. Action and Hamiltonian 


Starting from the Einstein-Hilbert action 

S[gf^,^]= J 

— yj det g—det g 


where is the spacetime metric with determinant g and R its Ricci scalar, the action can be 
written in Hamiltonian form as 

S[gij,TT"^]= j dx° J d'’'x gij - H) , 

where the tt*-^ are the conjugate momenta to the spatial metric components, defined by 


and H is the Hamiltonian density given by 


H = n^^gij - L = n^^TTij 




- gR- 


(3.1) 


This derivation is done e.g. in Appendix E of Wald |32) . 

The Hamiltonian density can now be written in terms of the Iwasawa variables and their 
conjugate momenta tt^, corresponding to /3“, and P\, corresponding to A/”® (note = 0 for 
a >i) which are defined as 


TTa 


dC 

5/3“ 


and P\ 


dC 

dN°'^ 


This gives 


H = ^ ca{N, P, 5,/3,5^/3,5,AA, , 

A 


(3.2) 


where G“^ = {S^^d - 1) - l)/(d -l),Af = (A/'“ ) and P = {Pi). The d x d matrix G“'’ is the 
inverse of Gab = - ^l^b^ which will appear later. In d = 3 dimensions they are explicitly 

given by 

/ 0 -1 -P 
{Gab) = -1 0-1 

-1 0 ; 


and (G“'') = - -1 


1 -1 -r 
1 -1 


,-i -1 


1 ; 


The sum in the second term of (3.2) contains the potential “walls” which will be discussed in 


detail in the next section. The derivation of (3.2), including the individual terms in the second 
part, is given in [Appendix A| 

The kinetic term JC only contains the conjugate momenta of the diagonal /3“ variables, the 
ones for the A/'“ are included in the “potential” term V. This makes sense because asymptotically 
the A/'“ tend to constants while the /3“ show linear behaviour, as will be demonstrated later. 
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4. The potential walls 


The structure of the potential term V in the Hamiltonian density (3.2) is crucial for the 
asymptotic behaviour. It is that of a sum, with each term consisting of a prefactor which, 
importantly, does not depend on ,5“ and an exponential term of the form exp(—2ajA(/3)) where 
loa is some linear form depending on the wall in question. Depending on the kind of wall, the 
index A can be a single or a multi-index 

The walls are split into two categories: The so-called “dominant” and “subdominant” walls. 
The dominant ones are defined as the minimal set of walls such that if their linear forms are 
positive, all the others are as well. Crucially for the billiard picture, the coefficients ca are 
positive for the dominant walls. 


The form of (3.2) allows the following “billiard” interpretation of the asymptotic dynamics 
(e.g. m) : A “particle” with coordinates moves through a Lorentzian space with metric Gab 
(from the kinetic part K. of the Hamiltonian) in a potential of the form V. The behaviour of the 
summands in the potential is dominated by the exponential terms exp(— 2a;yi(/3)). The /3“ can 
be decomposed as /3“ = p 7 “ with Gabl’^l'^ = — 1 and a heuristic argument, in analogy to the 
exact Kasner solution, gives p —>■ oo as r —>■ oo. In the limit, the potential walls become infinitely 
sharp as —2uja{I3) = —2pu}A{l) ±oo. As long as oja{I3) > 0 the potential is negligible and 
the 13°' evolve linearly. At the points where uja{I3) becomes positive the potential diverges and, 
because ca > 0 for the dominant walls, the particle is reflected. The subdominant walls do not 
influence the behaviour as they lie behind the dominant ones. 

This picture depends on the assumption that p —)■ cx) as r —> oo and that none of the 
walls vanish (either completely or asymptotically). The following does not depend on these 
assumptions, as the billiard picture will not be used. 

In the vacuum case there are two types of potential walls: The “symmetry walls”, coming 
from the kinetic terms of the off-diagonal metric components and the “gravitational walls” coming 
from the curvature term in the Hamiltonian density. The derivation of their exact form is given 
in [Appendix A[ here only the result is stated. 

4--1. Symmetry walls 


These come from the parts of the first two terms in the Hamiltonian density (3.11 which are 


not contained in the kinetic term K, in (3.2). The part of V containing the symmetry walls is 

(4.1) 




2g-2(/3'’-/3“) 


a<& 


where the multi-index A from (3.2) is {a,b) and runs over all a,b G {!,..., d}, a < b. The 
coefficients {ca) = {cab) are given by )^/2 and the linear forms {oja) = (wsym ab) by 


C^sym ab{/3) =/?'’- /3“ 


(4.2) 


The walls with the forms Wgym aa-i-i are the dominant ones among the symmetry walls, because 
if they are positive then > (3° Va and therefore all the other Wsym abiP), a < b are positive 
as well. 


4 . 2 . Gravitational walls 

These come from the curvature term in the Hamiltonian density ([ 33 . The gravitational 
walls split into two classes: 
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The contribution to V coming from the first class is given by 


E with a„,, = 2r+ E /3^ (4-3) 

a^b^c^a e^a,b,c 


i.e. the index A = (a, &, c) is a multi-index running over all a,b,c G d}, a ^ b ^ c ^ a. 

For d = 3 the sum in the expression for aabc{P) vanishes: There are only three possible values 
for the indices which are all occupied by a, b, and c, leaving no possible value for e. In this case 
only aabciP) = 2/3“ remains. 

The second class of gravitational walls has a more complicated form, their contribution is 
given by 

-E^^a(a2/3,a,/3,a,C',C)e-2'^“(« with /r„(/3)=E/3^ (4-4) 

a c^a 

and (all sums explicitly indicated) 


F,=- 2(/3“ - 2/3“ ,, 

+ E ( - ^(C'abr - 4Cl/3“ + 4/3;>,^“ - 

b ^ 

+ E [c'baC‘^ac - ” C\,C\J2 - , 


where the comma denotes the Iwasawa frame derivative Ca, defined in (2.2) and given in terms 
of partial derivatives as X^a = ■ Here H is a single index a G {1,..., d}. This term 

contains second derivatives of /3“ and A/'“ (it contains first derivatives of which contain first 
derivatives of A/’“ ) as expected from a curvature expression. 

The linear forms fia of the second class can be written as a linear combination of the ones of 
the first class, aabc, by 

dc — (o^abc 4 “ bXbca^ (2 . 


This means the first class of walls is dominant and the second subdominant. This is fortunate 
as the coefficients of the second class of walls, can be negative while those of the first class, 
(CV)^ /4, are always positive. 


1^.3. Complete Hamiltonian in Iwasawa variables 

The complete Hamiltonian density in Iwasawa form (equation (3.2) with the expressions for 
the walls inserted) is 


a<b 


4' 

a^b^c^a 


- E [ - 2(/3“ - 2/3“ + E ( - - 4G^,„/3“ + 4/3,",^“ - 

n • ' 


(4.6) 


- 2C\,I3'^, + 2/3“+ 2G“,., 


+ E ( C\aC\c - - C\,C\J2 - 2C\,P% 
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with CV = E^,k 2A/'“fc and G“'’ = {S^\d - 1) - l)/{d - 1) and where the 

derivative operator ‘\a” is defined as _q = 


5. Equations of motion and constraints 

For a Hamiltonian density of the form Hlq^Xjt), p(x,t), dxQ, d^q] the evolution equations are 
given by 


q{x,t) 

p{x,t) 


m 

dp ’ 


dq ™ ^^didmdnq) ' 


Here the variation is t aken after cho osing lapse and shift, which depend on the metric (the 
lapse is given by Vdet g). Appendix B shows that this does not change the resulting equations. 
In the case of the Iwasawa variable Hamiltonian (3.2) this leads to 


^r'^a — ^ ^ 


-"•"'(-TO'""™) 




dP\ 


drK =E 

A L 


_ J^y^p-2wAil3) ,a.f _J^yA_ 2 wa{i3)\ 

dUy ) 

^ ^\d{d,dkM\) J 


(5.1) 


where the components (uJA)a of the linear form oja appearing in the second equation are defined 
as (iXA)a = duJAil3)/dl3°-. 

The Hamiltonian and momentum constraints are 


where da 


H = 0, 

-\Ha ■■= db*\ + + C'^ac'^’^d - 


= 0 , 


{N-^)\d, and 


1 

p 


if 5 = a , 

2 -'' i- a - if 6 > a, 

|e-2(/3“-/3'’)Af“ P*j, , if6<a. 


(5.2) 

(5.3) 


(5.4) 


The Iwasawa variable form (5.3) of the momentum constraints is derived in Appendix C.l 


9 



















6. Fuchs theorem 


The following definition and theorem are by Choquet-Bruhat [33J Appendix V, p. 636], 
generalising the result of Kichenassamy and Rendall [53] . 

Definition 1 (Fuchsian System). A system of partial differential first order equations on V = 
M X K, M an analytic manifold which can be extended to a complex analytic manifold M, 

tdtu + A{x)u = tf{t,x,u, Dj;u), (6.1) 


with / linear in the first order spatial covariant derivative D^u, A and / extendable to holomor- 
phic maps in x and u (on M) and continuous in t G [0, T] is called Fuchsian if there exist a < 1 
and E > 0 such that := exp(A(z) logcr) satisfies 


sup 

zGM 




cr“ < E 


for te[o,r]. 


Lemma 6.1. A system of the form ^6.1\ ), with M, f as before, is Fuchsian if A is uniformly 
bounded on M x [0,r] with the real part of all its eigenvalues greater than —1. 

Theorem 6.2 (Fuchs theorem). A Fuchsian system has a unique solution u, analytic in x G M, 
in t and such that u = 0 for t = 0 in a neighbourhood of M x {0}. 


Replacing t in ( |6.1| ) by t' = gives 


p. ^t’dt'U + A{x)u = t'^ f(t{t'), X, u, Dxu ), 


and, as the eigenvalues of p ^A are simply those of A divided by p, the following corollary holds. 
Corollary 6.3. The theorem holds for 


tdtu + A{x)u = t^f{t, X, u, D^u) 
if all eigenvalues X of A fulfil Re{X) > —p. 

The condition that / be linear in the spatial derivatives D^u can be relaxed to admit an 
arbitrary analytic dependence by adding v := D^u as a new variable. Differentiating (6.11 gives 
an evolution equation for v, 


d f d f d f 

tdtv + Dj;A{x)u + A{x)v = tDxf = ^ ^ 


dv 


which is linear in Dj;V. Together with (6.1) this is a system of the form 

tdtU + A{x) = f(t, X, u, D,cu) 

for u = {u, v) and with A the block lower triangular matrix 


\D,,A a) ■ 


( 6 . 2 ) 


(6.3) 


The eigenvalues A of A fulfil 


0 = det(A — lA) = det(A — lA)^ 
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and are therefore exactly the eigenvalues of A. Therefore the system (6.3) fulfils the conditions 
of definition and is Fuchsian, provided / depends analytically on u, i.e. / depends analytically 
on DxU, and D^A is uniformly bounded. 


Corollary 6.4. The existence theorem \6.S\ holds for f depending analytically on D^u if D^A is 
uniformly bounded on M x [0,T]. 


A change of variables t = e J-rinlG.ll, with M a domain in M", gives the form of the 
theorem used here. 


Corollary 6.5. A system of the form 

drU — A{x)u = x, u, D^u ), 


(6.4) 


with A analytic in x and uniformly bounded, p > 0, f analytic in x, u and D^^u, continuous in t 
and bounded in t for t —>■ oo and with all eigenvalues X of A fulfilling Re{X) > —p has a unique 
solution u{x,t) with u{x,t) —?> 0 as t ^ oo. 

The matrices we will consider in the following will be constant and therefore the relevant 
conditions will be the boundedness of / as t —>■ oo and the condition A > —p on the eigenvalues 
of A. 

In order to obtain a more precise description of the decay of the solution, we dehne u = e'^'^u, 
0 < v < p. (6.4) becomes 


drU — {A + nl)u = e f{T,x,u{u),dxu{u)), 

which is again Fuchsian, as the conditions on / are unaffected and the eigenvalues of A+ i/l are 
shifted up to compensate the change in p. Therefore u = e'^'^u j 0 ^ 


= 0(e , V 0 < ^ < ^ . 


(6.5) 


7. Asymptotic evolution equations and differences 

7.1. Strategy and evolution equations 

We will use the strategy introduced by Kichenassamy and Rendall in [23] and used in m 
to prove the existence of solutions of the Einstein equations with non-chaotic asymptotics. As a 
first step we consider a simplified system of evolution equations, which is supposed to model the 
asymptotic behaviour, and which can be easily solved. Then we write down the equations for 
the differences between solutions of this system and those of the full one, following from the full 


evolution equations. If this system can be shown to be Fuchsian (i.e. if it is of the form (6.4)) 


then, by the Fuchs theorem, a unique asymptotically vanishing solution exists. This implies that 
a solution of the full system of equations exists, which asymptotically approaches the specihed 
solution of the simplified system. 

The constraints will be treated separately, in sections |7.2| and |7.3| 

Quantities relating to the asymptotic system will be marked with a subscript [ 0 ]; e.g. /Sf),,. 
The Hamiltonian of the asymptotic system is obtained by discarding all wall terms in the full 


Hamiltonian (3.2), leaving only 


1 


H[o] — ^G“^7r[o]a7r[Q]b. 


(7.1) 
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This gives the asymptotic evolution equations 


^r/3fo] = 2G“>]6, 
A/jo] i = 0 ; 


^T'^[0]a — 0 1 


"^[Oja — ^^abPt 

p i _ pi 

■^[0] a a ■ 


b 

o 1 


(7.2) 


with solutions 

/Sfo] =P> + /3o 
A/'[o]“ =-^0“ . 

Now, the differences ,5“, tfa, and P\ are defined as the real solutions min us th e asymp¬ 
totic ones (e.g., ,3“ = /3“ — Inserting them into the full evolution equations (5.11 gives the 


following equations for the differences: 

lgah-^^0, 


drTTa = ^ 


-ha, ( -2wA{0ioi)e-^^A{i3)'\ 

\d{diP‘>') J 


\d{didjl3°-) 


drN°‘- = ^ _^^2Lg-2i"A(/3[o])p-2u'A(^) 


dPi 


frP\=Y. 




dN° 

- djdk 


( dcA -2wA(0in]) -2 wa{P) 


dcA 


d{djdkJ\f°‘i 


^-2wAiP[0]) -2 wa(P) 


(7.3a) 


(7.3b) 


(7.3c) 


(7.3d) 


This system of equations is not directly in Fuchsian form, as the right-hand side contains second 

these can 


order spatial derivatives of the variables. By defining Bj := 5^/3“ and N° 
be expressed as first order derivatives. This is only possible if the j3 and Af equations (|7.3aP 


:= djAf^, 


and (7.3c I do not contain spatial derivatives, as otherwise new second derivative terms would 


appear in the evolution equations for the new variables. The (3 equation (7.3a) obviously doesn’t 
contain spatial derivatives while for the M equation (7.3c) the sum over the walls only includes 
the symmetry walls with coefficients )^/2, as the others are independent of P. 


The additional evolution equations for the new variables are given by 


1 


drB^ = -G“'’9,7fb, 


=Y,dj . 


(7.4) 

(7.5) 


To ensure that the right-hand side of the equation for decays appropriately we replace tt by 
TT defined as tTq := with e > 0. The evolution equation for Bj then becomes 


drB] 



(7.6) 
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The evolution equation for tTo, which replaces equation (7.3b) is then 


drTTa - ^ 


\0{0iP^) ) 

— f) d ( „-'^uiA(0m)„-‘2 wa(P)\ 


(7.7) 


with the additional term on the left-hand side and the exponential factor on the right-hand side 
coming from 

SrTTa = dr {e^'^TTa) = e^a + e'^'^dr^a ■ 


The full system of equations is now given by the 2d + d{d — 1) + d"^ + d‘^{d — 1)/2 equations 
([7^, (j^, ( [73^ , ([7^, ^ and ( f^ . 

The asymptotic behaviour of the terms on the right-hand side is dominated by the exponential 
terms exp(— 2a;^(/3[o])). If ojApQ]) is strictly increasing with r for all A, i.e. if uja{Po) > 0 for all 
A (for all walls) the system fulfils the decay condition required by the Fuchs theorem (Corollary 
6.5). Equation ( [F7t includes the exponentially growing term but as e can be chosen arbitrarily 


small, and therefore smaller than the minimum of uja{po), this does not affect the conditions. 

In order to be a Fuchsian system, the condition on the matrix A also has to be fulfilled. For 
this system the matrix A is given by 


/ Orf 

(G“^)/2 

0 d,d* 

0 d,d* 

0 d,d2 

0 d,d** 

Od 

eld 

0 d,d* 

0 d,d* 

0 d,d2 

0 d,d** 


0 d*,d 

Od. 

Od* 

0 d.,d2 

0 d*,d** 


0 d*,d 

Od. 

Od* 

0 d*.d2 

0 d*,d** 

0 d2,d 

0 d2.d 

0 d2.d* 

0 d2,d* 

0 d2 

0 d2,d.. 


^dtt,d 

0 d**,d* 

0 d**,d* 

0 d.*,d2 

Od.* / 


where d^, = d{d — l)/2, d** = dd* = d^(d — l)/2, 0^ is the d x d zero matrix, is the d x d* 
zero matrix and is the d x d identity matrix. The eigenvalues of this matrix are 0 and e. The 
condition therefore requires Re(0) = 0 > —uja{Po) and is fulfilled if uja{Po) > 0 '^A. 

In addition to the conditions uja{Po) > 0, the asymptotic Hamiltonian constraint, defined as 


Ho 


GabPtpl = 0 , 


(7.8) 


also constrains the values of the For vacuum in dimension d < 10 the conditions uia{Po) > 
0 VH cannot be satisfied together with the asymptotic Hamiltonian constraint Ho = 0 PO] . 
Therefore it is expected (e.g. m) that the generic solution in the vacuum case is chaotic. 

In d = 3 dimensions it is easy to see why the conditions are not compatible: The linear forms 
of the dominant walls (these are the only relevant ones) are 


Wsym 21 (Po) = pI - pI , UJsym 32(Po ) = pl - pl , Q!l23(Po ) = ‘2^pl 


(7.9) 


(two symmetry walls from (4.2) and one dominant gravitational wall from ( |4.3[ )). 
The asymptotic Hamiltonian constraint is 


Ho = -pIpI - PoPt - PtPo = 0 • 


1^3 


^2^3 


(7.10) 


The condition that the three linear forms (7.9) are greater than 0 implies Po>Po>pl>0 and 
therefore Ho < 0 which conflicts with the Hamiltonian constraint (7.10). 
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1.2. Asymptotic momentum constraints 

The asymptotic momentum constraints are obtained from the full momentum constraints 


lower triangular part ^^a[-] ^ diagonal part = —7^1,12, and discarding This gives 


(5.3) by splitting tt" 


(dehned in (5.4)) into a strictly upper triangular part a strictly 


1 

2 


Ha[ 0 \ —5b^[0]^a[-] “ 9^a’’'[0]o 


“ 2^l0fca'^l 


[0]a 


- -c 


+ ^[P\'cb'^[Cl]a[-] + C[0]^ac'^[o{d[-] 
[0]‘^ad’’'[0]d - 2 (^“^[0]'^)'^[0]d = 0 , 


(7.11) 


where da = (A/), ^Ya^i and are the structure functions of the asymptotic Iwasawa 

coframe, defined as in (2.3) but with AC® replaced with A/jpj “ = A/),“ . 

The only time dependent term in (7.11) is —da = -T{daPi)GdaPl- This term 

vanishes if GatPoPo = 0; i-®- T the asymptotic Hamiltonian constraint is fulfilled, because the 
matrix Gab is symmetric and constant: 


1 


{daPa)plGba = GbcT^ [ {daPa)pl + {daPl)Po ] = ^{9a{paPlGbe)) = 0 . 

=PoPo6=0 


1 


The asymptotic constraints are therefore preserved under the asymptotic evolution given by 

([TT]). 


1.3. Relationship between asymptotic and full constraints 

We want to show that if the solution (7.2) of the asymptotic evolution system fulfils the 
asymptotic constraints, the corresponding solution of the full evolution equations fulfils the full 
constraints (5.2), (5.3). 


The evolution equations for the full constraints coming from the full evolution equations (5.1), 
in Iwasawa variables, are 



drHa = Vai7 H- daQ , 

9 




(7.12) 

(7.13) 


with Ha (derivation in Appendix D). The right-hand side of (7.12) can be rewritten 

as 


E' 


,-2Ma(/3) 


daHa + 2{danHa ” 2 ( 9 ^ ^ d‘^)Ha 


with Paid) = J2b^ad’’ the subdominant gravitational wall forms. Defining H = e^^H, with 
77 > 0 , gives the system 


drH — T]H = ^ I 


,r;-2/i„(/3) 


daHa + 2{danHa - ^da^^ 


£7, 


drHa = e-"" ( VaH + -daP 


(7.14) 
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which is Fuchsian if 77 < 2fj,a(P). The term Vai^ in the second equation is equal to g da{H/g) + 
HO{C\^) = gda{H/g) + H 0{1) where the first part comes from the density character of H 
and the second from the connection coefficients in a non-coordinate basis. The system (7.14) is 
homogeneous and therefore the unique solution such that 




and Ha 


4 0 


(7.15) 


guaranteed by the Fuchs theorem is Ha = H = H = 0. We therefore need to check that (7.15) 
holds, i.e. that the constraints are asymptotically fulfilled. 

The differences between the asymptotic constraints and the full ones consist only of terms 
which vanish asymptotically: The Asymptotic Hamiltonian (El is exactly the part JC of the full 
Hamiltonian (3.2) which does not contain the exponential wall terms exp(—2 wq(/ 3)), which go 
to zero if the conditions u!a{Po) > 0 are fulfilled. The asymptotic momentum constraints were 
obtained from the full momentum constraints by discarding ff^^r,, = P\ /2, a > &, 


which is an exponentially decreasing term if the symmetry wall conditions (4.2) are fulfilled. 


This means that if the asymptotic constraints are fulfilled, the full constraints H and Ha vanish 
asymptotically. To verify ( |7.15 1 we still need to make sure that the definition of H does not 
change the asymptotic behaviour. This is the case as 77 can be chosen arbitrarily small, and 
therefore smaller than 2(/3“ — /3^), a > b, while still preserving the Fuchsian form of (7.14). 


This means, provided the solution of the asymptotic evolution equations satisfies the asymp¬ 
totic constraints, ( 7.15| ) is fulfilled. As the evolution equations for the full constraints are a 
homogeneous Fuchsian system, the unique solution which vanishes asymptotically is the zero 
solution. Therefore it suffices to impose the asymptotic constraints at one time (as they are 
preserved by the asymptotic evolution), to guarantee that the corresponding unique solution of 
the full evolution equations satisfies the full constraints at all times. 


8. Construction of the new class of solutions 


While generic solutions in the vacuum case are expected to be chaotic, there exist examples 
of vacuum spacetimes which show non-chaotic behaviour. As described in the introduction, 
all previous examples were at least 17(1) symmetric. These were constructed by starting with 
a symmetric ansatz for the metric, postulating asymptotic behaviour for its components and 
proving, via some sort of Fuchs theorem, that solutions with this asymptotic behaviour exist. 

Here, no symmetries of the metric will be assumed. The idea is to choose an ansatz for the 
AC® such that some of the walls in (7.3) asymptotically vanish. This means that their linear forms 


can be negative but the resulting exponentially increasing term is countered by an exponential 
decrease of the coefficients ca- 


8 .1. Ansatz and evolution equations 

The following ansatz is chosen for A7“ : 

AA“ (x4r) = A7„“ + e-^"A7,“ [x^r). (8.1) 

A/;“ is a constant, upper triangular matrix, with ones on the diagonal, which depends neither 
on space nor time. This ansatz for AC^ will cause the dominant gravitational walls to vanish 
asymptotically. It can be simplified to A/), = 1 by the space coordinate transformation a;* —> 
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y^{x^) defined by 


=x^ + Tv; ^2 + A/; ^3 x^ , 

y^ = x^+AfJ,x\ 


which does not affect the f3°‘. _ 

, /3“ and Ha are decomposed as before in (7.21, giving d* + 2d = d(d + 3)/2 functions of 
space Po\, pi and /3^ 

As before, the second derivatives on the right-hand side of the evolution equations are elim¬ 
inated by defining := and := and if is replaced by Ha := e^^ifa- The 

evolution equations for ,3“, Ha, A/’g® , , P“ and A^s“y are now 


9./3“- -G“'’7f6 = 0, 

dr^a - effa = ^ 




2uja(^[0]) -2u;a(/3) 


= E 


-f 

"■(£ 

- 

didj 

E 

dcA 

dP\ 

A 



dcA 


dU\ 

+ 


— 

djdk 


-2wa(P[o]) p-‘2wA{P) 


dcA -2wA(Pm) „-'2 wa{0)\ 


'-'3-''' i> 

dcA 


^- 2 wa{ 0 [o]) p- 2 wA{^) 


drN,% = 


4 \ a / 


drBi = e-^^-G^^djiTb. 


(8.2a) 


(8.2b) 


(8.2c) 


(8.2d) 


(8.2e) 

(8.2f) 


The additional term on the left-hand side of the A/), “ equation and the exponential factor on 
the right-hand side come from 


9. (e-^-A4“ ) = -7e-""A4“ + 
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The matrix A is now 


/ Od 

(G“^)/2 

0 d,d* 

Od,d* 

0 d,d 2 


Od 

eld 

0 d,d* 

Od,d* 

0 d,d 2 

0 <i,d** 


0 d*,d 

7l<i. 

Od* 

0 d*.d 2 

Od* ,d** 


0 d*,d 

Od* 

Od* 

0 d..d 2 

Od* ,d** 

0 ^ 2 ,d 

0 d 2 ,d 

0 d 2 ,d. 

0 d 2 ,d. 

0 d 2 

Odbd*. 


^dtt,d 

Od**.<i* 

Od**,d* 


Od.* y 


with eigenvalues 0, e > 0 and 7 > 0 and therefore fulfils the conditions for a Fuchsian system for 
all allowed decay coefficients on the right-hand side. 

To show the system (8.2) is Fuchsian, each term on the right-hand side has to be shown to 
be exponentially decreasing. 

To simplify the argument the decay rates as r 00 of the coefficients ca and their derivatives, 
with the ansatz ( 8 . 1 ) for A/"® , are now listed: 

dc, 


csym = {n^^r/2 = o{i) 


3c 

^^sym 


= 0 , 


syra ^ ^ 

= 0 , 


dcs 


dP 

dc 


sym 


_ _ Ofg-T-'-') 

5(9 ,^ ’ ’ 

Cs.d.g. = 0 {t'^) , 

3Cs.d.g 


didjakAf^^) 

dAf\ 

d^d.g. 


5(5,-5fcAA“ ) 

^^s.d.g. 


= 0 , 


5A/'“ 


Q^d.g. 

dP\ 


= 0 ( 1 ), 


= 0 , 


(8.3) 


d{djN\) 


= 0 {t) , 


5AA“ 

d^s.d.g 

d{djdkM\) 


cx /3“A(, = 0 (re-^"), 

= 0 ( 1 ). 


5c, 


s.d.g. 


dP\ 


= 0 , 


Here Cgym stands for the symmetry wall coefficients {Pd^Af^j )^/2, Cd.g. for the dominant 
gravitational ones, (O ^^)^/! , and Cg.d.g. for the subdominant gravitational ones defined in (4.5). 

As the ,3“ equation (8.2a) has no terms on the right-hand side we consider first the tTq equation 
(8.2b). The key terms are the exponentials exp(—2cij,i(/3)) and the ca and their derivatives. The 
overall space derivatives in the second and third part are innocuous as they can only bring down 
polynomial expressions in r from the exponentials. 

Beginning with the (relevant terms of the) hrst part, we look at the three 

kinds of walls (symmetry, dominant gravitational and subdominant gravitational) separately: 


symmetry The coefficients do not decay (see (8.3)). The whole term decays exponentially only 
if 

a^syni abiPo) P ^ )> 5, (^-4) 

i.e., as e is arbitrarily small, if the symmetry wall conditions are fulfilled (this means that 
Po < Po < • • ■ < pi)- 


dom. grav. The coefficients decay as e The whole term shows exponential decay if 


- 27 - 2aabc{Po) -I- e < 0 <t7 7 > -UabciPo) + e Va, 6 , c, 


(8.5) 


where the aabc are the linear forms of the dominant gravitational wall, defined in (4.3). 
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subdom. grav. The coefhcients do not decay, the whole term only decays if 


Aia(Po) > e Va, 


( 8 . 6 ) 


with fj,a the linear forms of the subdominant gravitational walls, defined in (4.4|. 


The second and third parts of the tt^ equation ( 8.2b[ ) include only the subdominant gravi¬ 
tational walls, as these are the only ones containing derivatives of /3. To guarantee exponential 


decay in these terms, the subdominant wall conditions HaiPo) > e have to be fulfilled, as in ( 8 . 6 ). 
The Ng “ equation, (8.2c), contains only the symmetry wall term, as the derivative of the 


other coefficients vanishes. Because of the exponentially increasing term exp( 7 T), decay requires 
that 


2a;. 


sym 


ab{Po) >7 ya> b. 


(8.7) 


The first part of the equation (8.2d I (containing A/”® derivatives of the coefficients) decays 
exponentially if the following conditions are satisfied for the different wall types: 


symmetry The symmetry wall conditions Wgym ab{Po) > 0 have to be fulfilled (as in (8.4|). 


dom. grav. 7 > —aabcipo) Va, 6 , c (as in (8.5)). 


subdom. grav. PaiPo) >0 Va (as in ( 8 . 6 )). 


The second and third part of the equation involve only the gravitational walls, as the sym¬ 
metry wall coefficients don’t contain any spatial derivatives of A/”® . Both the dominant and 
subdominant gravitational wall coefficients include first order (spatial) derivatives of A/”® , but 
only the subdominant walls contain second order derivatives. The dominant gravitational wall 
term gives the condition 

'Y >-2aabc{Po) Va, 6 , c ( 8 . 8 ) 

(2 because the decay of the derivative of the coefficient is only of order 0 (exp(—yr)) instead of 
0(exp(— 27 r))). The subdominant wall terms require Paijpo) > 0, as before. 

The A^s“j equation (8.2e) requires the same conditions as the A/),“ equation ( |8.2c ), as the 


overall space derivative only adds polynomial terms. 

Decay of the right-hand side of the P“ equation (8.2f) requires only e > 0. 
Summarising, the conditions are 


C^sym a&(Po) > e, 7 + «afcc(Po) >£, /Xa(Po)>e, 

2Wsym a&(Po) ^ 7 5 7 (Po ) 7* 0 , 6 7 0, 

for all indices a, 6 , c and the asymptotic Hamiltonian constraint 

Y.PoP'o=Q- ( 8 . 10 ) 

a^b 

Additionally, the asymptotic momentum constraint equation has to be fulfilled. This will be 
discussed in detail later. 
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8.2. d = 3 case 


In 3 + 1 dimensions the conditions (8.9) are explicitly (without removing redundant ones) 


^symiPo) ^ C 


pI-pI> 

Po - Po > e , 

P A Po > e > 

(8.11a) 

7 + CtabciPo) > e 


7 + 2pJ > e , 

7 + 2p2 > e , 

7 + 2p3 > e, 

(8.11b) 

Pa(Po) > e 


pI+pI> e, 

Po + Po > e , 

Po + Po > e, 

(8.11c) 

2‘tUsymiPo) ^ 7 


2(Po - Po) >7, 

2(Po - Po) > 7 , 

2(pApo) >7. 

(8.lid) 

7 + 2aabc{Po) > 0 


7 + 4pJ > 0 , 

O 

A 

(M O 

+ 

7 + 4p3 > 0 , 

(8.11e) 



e > 0 , 



(8.11f) 


and the asymptotic Hamiltonian constraint is 

pIpI + pIpI + pIpI = 0 . 

Let us show that they can be fulfilled simultaneously: 


( 8 . 12 ) 


The Hamiltonian constraint (8.12) gives 


1 PoPo 

Po=- 


pI 


pI 


(8.13) 


The co nditions (8.11a| follow from ( 8 .lid) by choosing e < 7 / 2 . Likewise, (8.11b) follow s from 
( 8 .lie) by choosing e < —2p\ (as 7 > 0 and p\ < 0). T he sec ond and third condition in (8.11c) 
follow from the first and from p^ > p^ > Q > p\ (from (S.llal), as does the first when inserting 
(8.13) and choosing e sufficiently small. The second condition in ( 8 .lid) follows from the first. 


The remaining ones are now 

m.m.{2{pl-pl),2{pl-pl)] > 7 > , 

(the last two conditions in ( | 8 . 11 e[ ) follow from > 0 and 7 > 0 ). 2 {pI - pi) > -dpj 


follows frompo > 0 after inserting (8.13). The last remaining condition, 2{pl—pl) > —4po, gives 

{V2- l)pl > pI. 

Summarizing, the conditions on p“ (at each spatial point) are 


Po > 0 , 0 < Po < {V 2 - l)p^ and pi =- 


pIpI 

Po +Pl 


(8.14) 


The remaining free parameters are the six functions of the space coordinates 


Po"a and 


8.3. Constraints 


In addition to the evolution equations, the constraints (5.2) and (5.3) also have to be fulfilled. 


The asymptotic Hamiltonian constraint was already included in the conditions discussed in the 
last section. We start by considering the asymptotic momentum constraints and then show that 
the full constraints are satisfied if the asymptotic ones are. 


Inserting the ansatz (8.1) (with = 1) into the asymptotic momentum constraint (7.11) 
gives 


1 


^ E - {GaoPla + PiaPlGdf) = 0 


Va, 


(8.15) 


b 

bl>a 


(all terms containing C'[q]“ 5 c vanish as the asymptotic limit of Af is constant in space). 
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This is equivalent to 


1^1,3 = -{pI+pD ^(Po.3+K.3 
Po\^ =‘^{G2cPl,2+Pi2PiGdf) 


'0.3 


+ Po,3 +/3o,3bo +Po) +/3 o,3(Po +Po)) 


Pr. 


1,3 


= -P. 


1,2 


2{G^,pl^^+Pt^^piGdf) 


(8.16) 


Given any functions j3l, P^\ and pi, pi, pi fulfilling (|8.14| one can determine Pi, P^\ and 


from (8.16), obtaining thus a solution of the asymptotic constraint equations. 


The full constraints are fulfilled if the asympotic ones are, following the arguments of section 


7.3 The condition that the (modified) evolution equations for the constraints (7.14) are of 


Fuchsian form requires that the subdominant gravitational wall conditions are satisfied, which is 
the case here (see (8.9)). To ensure the full momentum constraints converge to the asymptotic 
ones the symmetry wall conditions have to be fulfilled, which is also the case. Finally, the 


Hamiltonian converges to the asymptotic Hamiltonian if all terms after the first one in (4.6) 


vanish asymptotically. For the terms coming from the symmetry and subdominant gravitational 
walls this is ensured by the decay of the exponential terms, as for the general ansatz. For the 
dominant gravitational wall terms it follows from the decay of the coefficients, which contain 
spatial derivatives of A/"® , and the resulting inequalities ( 8 .lie). 


9. Analysis of the constructed solutions 


In the following we will need the asymptotic behaviour of the metric components of the 
constructed solutions. This is derived in [Appendix E 

From the results above there exist numbers 7 > 0 and ;/ > 0 so that 

/3“ = pIt + PI + 0(e-"") , = 51 + 0{e-G+Gr) , (gj) 

Setting 

o-po = Po + Po + Po > 0 , = Pi + Pi + Pi , (9.2) 

the leading order behaviour of the metric components, and those of its inverse, is 


900 

= - 

e 


+ 

O'') 


0 , 

9oi 

= 0 







Pii 

= e“ 

-2pIt- 


1 



912 

= e~ 

-2plre 


+ 

O'') 


0 , 

910 

= e~ 

-2p?rg 


+ 

O'') 


0 , 

923 

= e~ 

-2p^rg 


+ 

O'') 


0 , 


(I O''), 

5°* = 0 , 

g" =e2P°^+2^»(l + 0'^), 

5 I 2 = + O") ^ 0 

-13 ^ ^ 0 

-23 ^ _p2plr^2Pl (^2^ ^ 


where O'' := 0 {e~'''''), the are functions depending only on the spatial coordinates and with 
the behaviour above preserved under differentiation in the obvious way. The precise form of the 
K°‘i is given in [Appendix E 


9.1. Remaining coordinate freedom 

, with the choice A/), = 1 , constrains the 

remaining coordinate freedom. 


We wish, now, to analyse how the ansatz (8.1 
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In |12| it is asserted that transformations mixing time and space coordinates are prohibited 
by the choice of lapse and shift and the assumption that the singularity is approached as r —>■ oo. 
Presumably this should follow from the resulting equations 


dx"^ dx^ dr dr 


— det g 


[dry dx^dx^ 

\ Sr y df df 


= — det 


det g 


/ V r dx^ dx^ \ 
\ J dyl J 


(9.3) 

(9.4) 


(assuming a transformation t, x* —> f (t, x^), ?/*(t, x-^)). However, the assertion is not clear. 

An attempt was made to construct a Fuchsian system by starting from the transformation 
law of the Christoffels. Defining 


A3 := 


dx^ ’ 


and writing the transformation in terms of it gives 


dA° 


dx^ 




(9.5) 


We split A°‘p into A^-p = Ao°‘p+5A°'p with the assumptions Ao% = Ao“o = 0, Ao“& = Ao“f,(x®) 
and Ao°o = Ao°o(''")• Inserting this into the 7 = 0 equation of (9.5) gives 


drSA'^p = -5‘^5}drA,\ + {SA^rf‘^se + Ao%f\,){A,^p + SA^p) - (Ao“5 + SA‘^s)r^rg . (9.6) 


This system is unfortunately not of Fuchsian form, as many of the r“^.y diverge as r —> oo (see 
Appendix EI. 

Assuming nevertheless a transformation of the form 


(r,x*) (f(T,x*),y*(x'^)), 

and inserting into (9.3) gives dr/dy^ = 0 and therefore f = f(r). (9.4) then leads to 


9r \ 


’ 


which implies that the Jacobi determinant of the spatial transformation is constant (in space 
and time) and that f is an affine function of r. 

Starting with a general spatial coordinate transformation x* —^ y’’ {x^) the metric components 
gki in the new frame take the form 




(9.7) 


For gii this gives 


5ii = e 


1^) -w- 


_ P-2/31 


(9.8) 


with 0{X) = 0(exp(—2 pot)) denoting higher order terms. Assuming dx’^jdy^ ^ 0, which is 
necessary to preserve the conditions > pi > pi (8.11a), this implies — log(9x^/ dy^). 
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The equation for gi 2 is 




(9.9) 


Inserting exp (—2^^) from (9.8) yields 




(9.10) 


Requiring = Sf + 0(exp(—qr)) we obtain the condition dx^jdy'^ = 0. Similarly, from the 
gi 3 equation of (9.7), we get 


and therefore dx^/ dy^ = 0 . 

Using these conditions on the 322 equation becomes 


(9.11) 


-2/3' 

522 = e P 


fdx'^ 




+0{X) = 


=-2/3" 


(A/" 2 ) 


1 'l2^-2/3" 


(9.12) 


where 0{X) = 0{exp{—2plT)) denotes higher order terms and therefore /3^ = /3^ —log(9x^/ dy"^). 
The 523 equation is 


523 = 


dy'i QyS 


=0{X) 


(9.13) 


and yields 


= 


dx^ j dy^ 


■0{X) 


dx^ j dy"^ 

which implies dx"^ / dy^ = 0 . 

The conditions on the coordinate transformation are now 

x^ = x^{y^) , 
x^ = x^{y\y^), 
x^ = x%y\y\y^). 

Under such a coordinate change the asymptotic functions and /3“ transform as 


(9.14) 


(9.15) 


A/'o“ =K^=5t, 


^ 3x^ 

PZ=Pl, /3-=/3«_iog( — 


(9.16) 


As the remain unchanged, the conditions (8.14) are unaffected. Therefore 7 can be chosen to 
have the same value in the new coordinates, i.e. 7 = 7 . 


The transformations (9.16) reduce the possible isometries of the constructed solutions: Spatial 


isometries would have to be of the form (9.15) as otherwise the asymptotic evolutions would not 


match (a Mo 7 ^ 1 would not give an asymptotically diagonal metric and transformations which 
exchange the order of the /3“ would change the asymptotics of the diagonal terms). As pi and 
pI are only constrained by the inequalities 0 < < (v^ — l)po but otherwise free functions of 
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all space coordinates, which are not influenced by coordinate transformations of this form, we 
have the following proposition: 

Proposition 9.1. For a generic choice of the asymptotic functions and p^, the correspond¬ 
ing solutions have no (continuous or discrete) isometries f : M M of the form T{(j){q)) = 
T{q), \/q € M, i.e. involving only the spatial coordinates. 

9.2. Killing vectors 

In this section we wish to investigate continuous symmetries of the constructed solutions. 
This requires an analysis of the Killing equations: 

= X^d^g^^ + d^X'^g^, + d,X^g„^ . (9.17) 

We will seek Killing vectors of the form X = X'^d^ + X^di, with 

X" = 0(e-""), (9.18) 

and we will assume that the behaviour above is also preserved under differentiation. This ansatz 
is more general than the assumption of purely spatial isometries in the previous section, which 
would imply X'^ = 0. It does, however, only include isometries which can be described by Killing 
vectors, i.e. continuous but not discrete ones. 

The ri killing equation states 

.^xgri = drX^gji + diX^grr = 0 . 

Contracting with gives 

= -g^^grr^^X^ = 0{e-(‘^P° ~ 2^°+-)-) = 0(6"""), 

i.e. the r derivative of the spatial components of the Killing field decay exponentially. 
Considering the rr component of the Killing equation gives 

= X^d^grr + X'^drgrr + 2grTdrX'" 

= J + JW - d^xA (1 + O'') . (9-19) 

'-V- 

=0(e-''^) 

The term of order one inside the bracket is non-zero for large times unless 

Xld.Up^ = 0 and Xld.ap^ = 0 . (9.20) 

The ii component of the Killing equation gives 

■^x9ii =X'''^kg^^ + X'^^rg^^ + 2diX^gki = 

-2X„^9fc(p*r + /3:)- (9.21) 

-f 2a,ds::o(e-2(Po^+^°)), 

which is certainly non-zero unless the highest order term, containing p(t, vanishes. This requires 

X^dkPl = 0 . (9.22) 
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There are no solutions fulfilling these conditions in general: The second equation of (9.20) 
and (9.22) can be combined into = 0, with A containing the derivatives of pi, pi and . 

If the determinant of A is non-zero the only solution is = 0. If this holds in the neighbourhood 
of a point, the Killing vector vanishes everywhere. We conclude that our solutions will not have 
any Killing vectors of the form (9.18) in general. 


Proposition 9.2. For a generic choice of the asymptotic functions pi, pi and j3f the corre¬ 
sponding solutions contain no Killing vectors X = X'^dr + X'^di satisfying X'^ = for 

any e > 0. 


Considering now a general Killing vector field X^, satisfying 


V ^X,^ V— 0, 


(9.23) 


we start with the expression VrVaX^. Using the definition of the Riemann tensor we obtain 


XrXaXjS = gp^V aX'^ = XaXrX-^ -\- RpjraX°' . 


Using the antisymmetry of the Riemann tensor and the first Bianchi identity gives 

XrXaXp = VaXrXj — R-^praX'^ = VaXrX-y {R-yra/S + R'fa/3T)X'^ . 


Applying (9.23) and using the definition of the Riemann tensor again leads to 


R-ycfi-rX-' 

XrXaXp = R-yrapX'^ — X aX pXr —XpXrXa + X rX pXa 

= R-frapX'^ — XaXpXr -\- XpXaXr — XrXaXp 
= ‘^^R'fTa/sX^ — XrXaXp , 


and therefore 

XrXaXp = R-yrapX'^ . 


Introducing 

FaP = -^{XaXp — XpXa) = XaXp , 

we thus have the following system of equations for the pair {X, F) if A is a Killing vector: 


XrFap = Rara/sX'^ ■ 


(9.24) 


This has the general structure of a Fuchsian system, but does not fulfil the conditions given in 
section By redefining some of the A“ and F^p the equations can be brought to a form where 
all terms on the right-hand side decay exponentially, but the equation for X^ contains the term 
FqqA^ on the left-hand side, which gives a negative eigenvalue —{pi +pI+pI) in the matrix A. 
To still satisfy the conditions of the Fuchs theorem all terms would have to decay faster than 
exp(— t(po +pI -\-pl)), which is not possible. 


9.3. Relationship with previously known solutions 

The first class of vacuum spacetimes for which asymptotically simple behaviour was shown 
was the polarized Gowdy class Emu]. This is a class of solutions containing two commuting 
spacelike Killing vector fields (the polarization condition) with constant t hypersurfaces which 
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are compact without boundary and orientable (the Gowdy condition). The topology of spacelike 
slices of these spacetimes is constrained to one of x S^, or a Lens space L{p, q). In the 

following only the case will be considered. The metric for this case is given by 

ds^ = + d9^) + tie'^dx^ + e-^dy'^), (9.25) 


where a and W are functions of t and 0 which are 27r periodic in x. 

Redefining the t coordinate as t = e~'^ transforms the metric to 

ds2 = ^W-t^^ 2 ^-W-T^y2 ^ 

which is directly in the gauge used here: The shift vanishes and the lapse is equal to the 

determinant of the spatial part of the metric. Comparing with (10.11 shows 

=e2“^/3i = -a, 

e - =e^ 


„- 2 / 3 " _ ^ ^2 ^ T-W 




/?" = 


2 ’ 
r+ 1T 


M\ = N\ =M\=d. 


The results of Chrusciel, Isenberg and Moncrief show that solutions of this form are parametrised 
by two functions of the 6 coordinate, tt and w, appearing in the asymptotic expansion of a and 
W [51]. The expansion is 


W = 7r(T — To) +UJ + 0(Te ^’’), 
a = (1 - 7 t^){t - to)/ 4 + a + 0(re“^'^), 

where tq is a constant and a a function of 0 which can be determined from tt and uj. This gives 
for the Pa 

I/O O 1/^ \ O 1 


pi = - 1): 


Po=o(l-7r), Po 


= 2 ^' 


(9.26) 


These satisfy the asymptotic Hamiltonian constraint ( |7.8[ ) but not necessarily the inequalities 
(8.14). These would imply ^2 — 1 < tt < 1 but there are solutions of the form (9.25) for any 


function tt. The assumption that the metric coefficients only depend on the 0 space coordinate 
and Af® = causes some of the potential walls to vanish: The coefficients of the dominant 
gravitational walls are proportional to A/”® ^ and therefore vanish identically. For A/”® to be 
constant, P\ has to vanish, which causes the symmetry walls to vanish. The coefficients of the 
subdominant gravitational walls (4.5) contain terms which are not proportional to A/'“ ^ . These 
do, however, contain spatial derivatives of the j3, most of which are zero here. As /3 only depends 
on t and 0, only one of the walls, with linear form = /3^ + /3^, remains. Therefore, assuming 

that the metric coefficients depend only upon 0 and that Af^ = 1, the only conditions left in 


(8.14) are 


Po = - 


pIpI 

pI +pI 


and Po + Po > 0 . 


(9.27) 


These conditions are satisfied by (|9.26[) which implies Po + Po = 1 and pj = —plpf 


(One should note that not every solution satisfying (9.27) is of polarized Gowdy type, as Po 
and Po can still be independently specified.) 

More general (non-Gowdy) symmetric spacetimes have also been shown to exhibit simple 
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asymptotic behaviour. These take the general form 

ds^ = + dx^) + e'^^{dy + Adz + (Gi + AGzjdxf + e-^^t^(dz + Gzdxf , 

with 77, [/, a, A, Gi and G 2 depending only on t and x [26]. To obtain simple behaviour either 
polarization, corresponding to A = const, or half-polarization, corresponding to a restriction on 
the asymptotic behaviour of A, has to be assumed. In both cases the resulting spacetimes are 
not contained in the class constructed here. The functions Gi and G2 tend to constant (in t) 
functions of x, but they appear in the Af® in the Iwasawa decomposition. This conflicts with 
the assumption Af —> 1 (or —>■ const) made in constructing the new class. In this sense the new 
class is therefore more restricted than the polarized and half-polarized classes. However it 
includes free functions depending on all space coordinates, not just one. 

The Killing vectors of the Gowdy and general T^-symmetric spacetimes are of the form 
considered in section [9^ as they do not include derivatives with respect to t. These are therefore 
in general not present in the class of solutions constructed here. 

10. Conclusion 

We have constructed a new class of four-dimensional (analytic) solutions to the vacuum 
Einstein equations which show asymptotically simple behaviour near a spacelike singularity, 
approached as t —> 00. The metric takes the form 

ds^ = -e-2^»^“dT2 + Af“ dxMx^, (lO.I) 


with /3“ and A/’“ depending on all coordinates t, x* and behaving asymptotically as 
/3“ = /?“ + tpI + 0(e-"") and Af“ = + 0(e-(^+")"), 


( 10 . 2 ) 


where 7 and v are positive constants. 

The class of solutions includes three completely free functions of all space coordinates, /3o, 
/3o! Po\ iPo\ does not appear in (lO.I) and (10.2) but influences the exponentially decaying 
terms, as detailed in Appendix E) and two functions, and p^, also depending on all space 
coordinates, which are constrained by the inequalities 

p3>0, 0<pl<iV2-l)pl. 

The Kretschmann scalar of the solutions behaves as 

K = = (G^ -b C>(e-'^^))e^4(P=+P°+P°) , 


with G^ a positive constant, dehned in 


E.12 


As pI + pI + Po > d the curvature tensor grows 


uniformly without bounds along all causal curves. 

For a generic choice of the free functions, the solutions have no continuous or discrete spatial 
isometries, and no continuous isometries described by Killing vectors satisfying X'^ = 

0(e-"^), e > 0. 

The construction is unaffected by the presence of a cosmological constant, as demonstrated 
in [Appendix F[ giving the same asymptotic behaviour for all values of A. 
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Appendix A. Derivation of Iwasawa variable Hamiltonian 


Here we will give the derivation of the Hamiltonian density in Iwasawa form (4.6), from the 
standard form of the Hamiltonian (3.1). 

In the following the spatial metric and its inverse will be used in their Iwasawa forms 


Appendix A.l. Kinetic and symmetry wall terms 

The kinetic term K. and the symmetry wall term come from the first two terms in the Hamil¬ 


tonian (3.1). These are 




(A.2) 


The conjugate momenta in Iwasawa variables, and , can be expressed in terms of tt*-^ as 

TT, = ^ = ^ = -2e-2/3“Af“ , (A.3) 

5/3“ * J > 


and 


P\ = 27r*^e-^'^“Af“ 


(A.4) 


We start by considering the first term in (A.21, n'^^TTij. Lowering an index in the first component 
and raising one in the second (using the Iwasawa form (O of the metric) gives 


The double sum can be split into a diagonal and off-diagonal part 

+ ■ (A.5) 

a a^b 


Raising the index k on in the first, diagonal, part leads to 
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where in the last step the definition of TTa, (A.3), was used. This, together with the second part 
of (A.2), gives 



a 


^ \ a ^ a ) 

^ V a ^ a / 


which is the kinetic part /C of the Hamiltonian (3.2). 

Raising the index k in in the second, off-diagonal, part of (A.5) gives 


b^a b^a 


This is symmetric in a and b and can be written as 


2^ AA';- 

a<ib 


a<b 


which is the potential term coming from the symmetry walls. 


Appendix A.2. Gravitational wall term 

The gravitational wall term comes from the term —gR in the Hamiltonian (3.11. 

We will calculate the curvature scalar in the Iwasawa frame. The Cartan formulas for the 
connection one-form uj‘^^ are 


d6l“ + ^ 0;“^ A 0'’ = 0 , 

(A.6) 

d'7afe — ^ab “t“ ^ba ? 

(A.7) 

where 

lab = 5ab exp(-2/3“) = 6abAl , 

with Aa := exp(—/3“), is the metric in the Iwasawa frame. We will also 

use the definition of the 

structure functions 

d6»“ = A 0G 

2 

(A.8) 

can be obtained by considering the expression (no summation) 

Alde^{ej,ea) + A|d6'-^(eb, Ca) - A^d6»“(ej, Cb). 

(A.9) 

Using (A.8| this is equal to 

A2f-fb 1 A2f-fa 

^6^ ja ba ^ jb • 

(A.IO) 
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Starting again from (A.9) but using ( |A. 6 [ ) gives 

- u;\{ej)) +A'^^{uj\{ea)-uj\{eb)) - AKuj'^jicb) - uj\{ej)) . 

Lowering the upper index on w with 'jab and using ( |A.7[ ) in the form ijJabi^c) = — Wba(ec) + 
i5ab(A^) c (with , c denoting the frame derivative by Cc) we obtain 


^jb{Aj)^a + ‘2uJab{e.j) 5ab{Ag) daj{A^)^b- 


(A.ll) 


Setting (A.10 1 equal to (A.ll) and raising one index using 7 “^ = gives a;“f, as 

‘^°‘bi^j) = + 6ab{A'^)j + 5aj{A^a) ,b “ 5jb{A'^)^a^ . (A.12) 

The curvature scalar can now be computed using 

n\ = duj\ + Y,bo\Au:^b 




ej 


(A.13) 


as 


= E ^ E A Cb) 


e,f,a,b 




(A.14) 


We start by calculating Aj, ^da;^(ea, Cb). This gives (with summation over all indices which 
occur more than once) 

^b^<i^^bi^a,eb) = 

^b"^ [ea(w“h(eb)) - ebiuj\{ea)) - C‘^^b^\{ec)] 


= A 


-2 


^ 2^4 (^^bC'^ab + 2Sab(Ab)^b — (Ab),a 


+ i‘^{^b) ,aC^ab + ‘^Sab{Al) ^b,b — {Al)^a,a + 2AlC^^b,a) 


^ 2^4 (2^aC'“ab + i^a) ,b) — 7^ ,bC°‘ab + i^a) ,b,b + 2AlC‘^^b,b] 


2 A2 


C 


+ ^cC'^afc + ^l^°'cb + 5ab{Al)^c + ^aci^a) ,bScb{Al) ^a) 


2 A 2 


= -C" 


ab 


{AD,, {Al),b{Al),b ^ {AD,,{AD,. 


A4 




‘^AtAl 


+ 2 C'*’ 


{AD,, ^ {AD,b,b 


^^AlAl 


At 


{AD,a,a , A,a , {AD,b{AD,b C^bO 


AlAl 


+ 2 


Al 


A 2 


2A4A2 


lb ac _ c _ ^ 

A2 W abJ 2A2A^ 


(AD,c 

bb 2^4 • 
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The second term, bc^b ^ e^), gives (again with all sums implied) 

^fr^(^“c(ea)w=b(ef,) - w‘^^{eb)u}'^biea)) 


4^6 L 


+ 2 C“ , + 


Al 


(Al),. 


Or^ ill _|_ 2 <^eb(^D,b _ 

cb ^2 42 42 

-^1, 


A? 


I /^b b I /^c c /oa I 
\ ^ ar. ^ ah ao ^ r.b ' 


Sac{Al),b , 5ab{Al),e Sa{Al),a 


ac A 2 ' a-b'A 2 


+ 


42 42 

I ^b ^b y^i 

cb ca ~ ^ 

AA-fA AAfA 


cb ' A 2 ' 42 42 

^cb(^^),a , S,Ml),b Sba{Al), 


= C“ C",,—+ 6 -“ 


ac cb j ^2 ' ^ a,b ^4 


' Al ' A 2 
(Ag),, , {Al),b{Al),b {ADAAD,. , 






42 

(^a. \2 ~^a 

^ 4412412 


+c“,c 


1 


(^g).b(^D,b , (^D.a(^D.a , (^D.c 


2416 


+ ■ 


^AlAt 


^a v~ ~f> / 1 ^; _ ^0 ,a 

^ c„ 42 42 ~ ^ ab 42 42 ■ 

-r±uA±^ 


Adding the two expressions and substituting {AA,b = —‘^A^I3‘\ we obtain the curvature scalar 


^ = -iE(c'“cbf 


AlAl 




A 


b , , 


a a L 

b / 3 'a 


_ 4(7“ , 1^+4 

cf> ^2 ^ ^ ^2 


_ n^b ' _ 

ab <2 42 

■^a -^a 

a /:?6 


A 2 


+ 2^+C“,,E (A.15) 


/ 1 ^ fi^ 

I \ '' [ __/^a ^6 _ ,cr',c _ ^ ^b _ ^ 

^ / A (^^26^ ac'^ cb j ^2 2 A^^ ca ^2 


Multiplying this with —g = — exp(—2 /3“) gives the gravitational wall terms (4.31 and 


(4.4) in the Iwasawa variable Hamiltonian. 


Appendix B. Iwasawa evolution equations and Einstein equations 


To obtain the evolution equations (5.1) the variation was taken after choosing lapse N and 
shift ]Va with the lapse given as -^/det g, i.e. dependent on the metric, and the shift vanishing. 
The general Hamiltonian density, with lapse and shift still free, is 


-A+(detg) Wij - - (det g) ^(7r))2 


iJ = \/det g\ N 

- 2A4V4(detc()-i/27ry)] +2V4(det5)-i/2A,-7r*^‘) 


(B.l) 


(Equation (E.2.32) in Wald [52]1. 

Taking the variation of the Hamiltonian H = J Hd^x with regards to N and Na gives the 
Hamiltonian and momentum constraints, respectively: 


£ TT 1 

— = -R+ {detg)-\^^Tr,j - -{detg)-^{TTlf, (B.2) 

£ pT 

— = V, ((det5)-i/2^^*) = V.ttA . (B.3) 
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Varying now with respect to and gij gives the Einstein equations in Hamiltonian form as 


9ij = = 2(det5) V + 2V(,V,) , (B.4) 

+ ^/V(det5)-1/"5*1 (B.5) 

- 2V(det 5 )-i/2 - ^7r'=fc7r*^^ - \/det g{\7^\7^N - g^^V^VkN) 

+ Vfc(V'=7rd) - 


(Equations (E.2.35) and (E.2.36) in Wald [52]1. 

Choosing Na = 0, either before or after varying, just removes the terms containing Na 


the evolution equations. Choosing N = -y^det g before varying adds an additional term in (B.5). 
This term is, however, proportional to —R + (det g)~^TT'‘^TTij — |(det g)~^{'n'l)‘^ which is zero, by 
the Hamiltonian constraint (B.2|. 

The terms in (B.5) which contain covariant derivatives of the lapse also vanish, as the deter¬ 
minant of the metric is covariantly constant. 

The transformation to Iwasawa variables is a point canonical transformation and therefore 
doesn’t change the equations. 


Appendix C. Derivation of Iwasawa variable momentum constraints 

In this section we will give the derivation of the momentum constraints in Iwasawa variables 
and the definition of their asymptotic equivalent, following section 3.2 of m- 


Appendix C.l. Full momentum constraints 

We start with the momentum constraints in the form 


= 0 


(C.l) 


(see e.g. equation (E.2.34) in Wald [H2]b 

The calculation is simpler when done in the Iwasawa frame (2.2), so we first calculate the 
Iwasawa frame components of in terms of the Iwasawa variable conjugate momenta and 
tTq. These are denoted by 7f“^ and defined as 

:= 7r(0“,0'’) = . 

Starting from gijir''^ = +-^“iPo writing the left side in Iwasawa variables gives 

= p'^TTa + ./V“ . (C.2) 


Using the diagonal form of the metric in the Iwasawa frame, jab = exp(—2/3“)Jat,, we obtain 


^ a ~ 


(no summation), 


(A/' 


— l\i 


= -P' 

c'‘ a 2 ^ 


(C.3) 

(C.4) 
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for i > a, 










by comparison. This can be written as a matrix equation 


(A/" ^)*c(+)''‘ a - 2 ^ a[-] a(+) 




+ X\ 


(C.5) 


(for all i and a) where the subscript (+)/(—) designates an upper/lower triangular matrix and 
[+]/[—] a strictly upper/lower triangular one. The matrix is defined by this equation. 

Multiplying (C.5) hy gives 


(+) ^ “ [-1 + (+) ^ “(+) ■ 


The strictly lower triangular part of j_j, is given explicitly by 

1 


^ a[-] = with 9{x) := 


0 if X < 0 , 
1 if X > 0 . 


Because of the symmetry of 7f“*' this also gives the upper triangular part via 

=0{a - = 0{a - 

=e-^iP^-P^)^\0{a -b) = . 

Finally, also including the diagonal term from (|C.3 ), we arrive at 


-TTb 

7f ^ i i A/'*’- P* 


for a = b, 
for 5 > a , 
e-2(/5“-/3'’)Af“ for a >6. 


(C.6) 


Writing the momentum constraint (C.l) in the Iwasawa frame and expanding the covariant 
derivative gives 

, (C.7) 

where the last term comes from the fact that is a tensor density of weight 1 and T^j^, are the 
connection coefficients in the Iwasawa frame, given by 

1 


^ &c 2 ^ ^ 0 i.9tTC,b ffi 9bt7,c 9bc,{7 P^tjbc “f P'bac T P'ccrb') 

(T 

= {dac{e~^^‘‘),b + ^ah(e“^^“),c - (5bc(e“^^*’),a 

- e-2/5“CV + , 

(no implicit summation) with the comma denoting the derivative in the Iwasawa frame. 
Inserting (C.8) into (C.7) gives 

- l/3\^TTb , 


(C.8) 


(C.9) 


which is the momentum constraint (5.3). 
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Appendix C.2. Asymptotic momentum constraints 

The asymptotic momentum constraints are obtained from the full ones by discardi ng th e 
contribution of the strictly upper triangular part of tt^^: Indeed, this part (the last line of (C.6)) 


contains exponential terms which vanish asymptotically if the symmetry wall conditions are 
fulfilled. 

Splitting into = ^^a[+] + ^^a[-] ~ <^6’’’a/2, discarding ^^a[+] inserting into (C.9) 
gives 




^ ^ ^^a[-\,b ~ 2^a,a + C'‘^cf)'^^a[-] + C'^ac^^&[-] 

- \c\a'^a - \c^abT^b - , 


which is the asymptotic momentum constraint (7.11). 


Inserting the ansatz ( |8.1[ ) (which implies = 0 asymptotically) and the asymptotic evo¬ 
lution of the /3“, = PoT -|- /3“ gives 

“ 2'^[^]a,a “ 2 ^ Aol^.a^fOlfa 

h b 

= E 1^0 ■ Po \b - E ^<^»Po,a - E(^°.a + rpljG,,P: ■ 


b>a 


b,c 


The last term, which contains a time dependence, is zero if the asymptotic Hamiltonian constraint 
GbcPoPo = 0 is fulfilled. Expressing the Iwasawa frame derivative in terms of the coordinate 
derivative ft leads to 


^ ^ E-^o(A/; -ybd,p^ - ^(Afo -y, | ^ Ga^dy,+ ^ G^fpioyt |, 

dj 


6 >a 

hj 


which is (8.15) with arbitrary (constant) A/, 


Appendix D. Evolution equations for the constraints 

Here we will give the derivation of the evolution equations for the constraints in our choice 
of gauge. Our treatment is similar to, but not identical with appendix A of |dl] . 

The hrst order action corresponding to the Einstein equations is given by 


TT*’, A, A*] = / dx^d‘^xigiy^ - NH - N^Hi) 


(D.l) 


with A the “rescaled lapse” defined as N = A/y’p and A® the shift vector. In our choice of gauge 


A = 1 and A® = 0. Erom (D.l) the equations of motions, the Hamiltonian constraints and the 


momentum constraints can be obtained by varying with respect to pij, A and A® respectively. 

We will compare the resulting equations with those coming from variation of the standard 
Einstein-Hilbert action Sh = / d^Xyt—gR^ which is (neglecting boundary terms) 


6Sh = - d^xV^G^^'Sg, 


(D.2) 
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with G denoting the Einstein tensor. 

As the spacetime metric is defined, in terms of N, Ni and gij, as 




(NkN'^-N^g N,\ 

\ g^J ’ 


the variation of Sh with respect to gij, N and Ni following from (D.2) is given by 


^ = Ng(-&^ + N^gG°°g^^) + 0{N^), 

^9ij 

^ = 2G°°N^g^ + 0{N^) , 

5N 

^ = -2NgG^° + 0{N>^), 


(D.3) 

(D.4) 

(D.5) 


where 0{N'^) denotes terms proportional to which vanish in our gauge. Here ^/—g = Ng 
and 5g = gg''^ 6gij were used. 

From the first order action (D.l) we obtain 


5N 


6S 

SN,. 


= -W . 


(D.6) 

(D.7) 


Identifying (D.4) with (D.61 and (D.5) with (D.7) yields 


H = -2g^N'^G°° + OiN'^) = -^Gqo + 0{N'^) , 

N'^ 

W = 2gNG^° + 0{N^) = -^g^^Goj + 0{N'^) . 


(D.8) 

(D.9) 


Using (D.8) to rewrite the equations of motion SSn/Sgij = 0 from (D.3) leads to 


= - — Hg^^+0{N'^). 
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(D.IO) 


We now consider the (vanishing) divergence of the Einstein tensor VyG''^ = 0. Using the 
identity 


— 

V=g ’ 


this can be rewritten as 


0 = = d,G\ + 

= d,G% + - do9ui.)G\ 


Expressing the components of the Einstein tensor using (D.8), (D.9) and (D.IO) gives for g = 0 
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(with dr = do) 


O(N^) = ^ ( d,{9ooHy2) - dr ( 500 


H 

' 2 ^ 




TTi /VI FT 

= - :^ddNg) - —d,W + —drH + -^drN -^ 

2Ng ^ 2 2g 2Ng Ag^N^ 

The last three terms cancel as dr g = gg^^drgij- 

Doing the same for the fj, = i equations, we obtain 


^ dr{N^g) + ^g^^drg.,. 


0{N^) = ^ - ^d, 
^ ^ 2Ng 2 * 


TT TT 

-^d,g - ^d,N. 
2g^ N 


(D.ll) 


(D.12) 


As is a scalar density of weight — 1, iJ a scalar density of weight 2 and Hi a tensor density 
of weight 1 their covariant derivatives are given by 

V,N = d^iNy^)/y^, V,H = gd,{H/g) and = d,H,. 


The divergence of can therefore be expressed as = di{g^^Hi). Inserting this into (D.ll) 
and (D.12) and applying our gauge choice = 0, iV = 1 yields 


drH = g^,W + Wd^g , 

drH, = V,H + -d,g . 
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(D.13) 

(D.14) 


The last term of equation (D.14) seems to have been overlooked in the derivation of m- It does 
not affect the arguments of section 7.3 concerning the relationship between asymptotic and full 
constraints. 


Equations (D.13) and (D.14) can be expressed in the Iwasawa frame as 

drH = 5Vai7“ + H'^dag , 

drHa = VaH + —dag, 

9 


(D.15) 

(D.16) 


using iJida;® = Ha9‘^ = HaN°‘i dx® with Ha the components of the momentum constraint in the 
Iwasawa frame and da = {JV~^)^a^i Iwasawa frame derivative. 


Appendix E. Asymptotic behaviour of the metric and curvature 


The ansatz (8.1) and the Fuchs theorem give for the asymptotic behaviour of the matrix JV 
A7“ = 6^ + = St + . 

An explicit first order term for the off-diagonal elements can be found from the evolution 


equation for , (8.2c). This equation is 


= yAA/, + ^ ^g-2^.4(/3[o])g-2^^(« ^ 


(E.l) 


The sum over the potential walls includes only the symmetry walls, as the coefficients of the 
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others do not depend on . Their coefficients are Cab = which gives, inserted 

into JETI), 


-2(/3‘=-/3“ 


(E.2) 


into (E.l), 

=7A4“ 

c>a 

Inserting the asymptotic behaviour of P and M gives 

+ e"" 5] {{PA+Kji^^ + )(5 + . (E.3) 

A ^ ' B " C ^ 

We consider first term A. As Po^a is of order 1 while P\ decays as by the Fuchs theorem, 
the first term dominates. For term B the situation is similar. The Sj term is of order one for 
c = j (the other terms in the j-sum are irrelevant as they share the same decay rate) while the 
second term decays as 
This leaves 

drJ^A = . (E.4) 

c,c>a ^ 

For C the situation is different, as the sum over c also influences the factor exp(—2(/?'^ — /3“)). 
We will consider the three components of “ separately. In the case A/), \ , i.e. a = I and i = 2, 
only the summand with c = 2 survives in the sum, as A4° = Afs% is zero for c = 2 and c = 3 
because of the upper triangular form of A/),. The evolution equation therefore becomes 

drAf, \ =-fAfs \ + 

Integrating this and using the fact that = 0{e~’''^) from the Fuchs theorem to eliminate 

the integration constant gives 

- Pq 1 e ^-r( 2 pl- 2 pl-'y} _|_ Q/g-'r(2p^-2pi-7+1/) \ 

2 (Po-Po) 

- 


K\=-'- 


(E.5) 


=:ifi 


For Afg ^3 the situation is similar: Afg ^3 vanishes, leaving only the Sf term and giving 

d^Mg% = -fMg% + e^^P,\ + 0(e-"(2p^Vo-7+^^)) ^ 

and therefore 


K%=- 


T~i 3 —2f/3^ —5^') 

2^ g-T(2Po-2Po-7) 

‘^{pI-pD 


(E.6) 


=-.kK 


For Ng ^3 two summands remain: For c = 2 the term with Ng\ and for c = 3 the one with <5°. 
The evolution equation becomes 

a.W ‘3 =7W‘, 
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Inserting A/'s ^3 from (E.6) gives 


drM, 


s 3 







Pq 1-Pq 2 \ -T(2p/!-2pi-7) I p-T{2pl-2pl-'^+v)\ 

2pl-2pl) + ^ ^ 


and after integrating 


K 


s 3 


= e- 2 (/ 5 o-/ 3 o)fp 3 _ Po\Po\ \^^ - _r{ 2 pl- 2 pl-^) ^(j(^^-T{ 2 pl- 2 pl-^+^)y 

V 2pi - 2pi J 2pi - 2pl 

= -.K\ 


by 


From equations (E.5), (E.6) and (E.7) the full matrix A/'“ and its inverse {Af are given 


= (K\ + 

Af\ = {K\ + 0‘')e-P'^P--'^P-'> , 

M\ = (X 23 + 0 ")e-"( 2 p^ 2 p^) , 


= -{K\ + 0")e-"(2p?-2Po) , 

= -{K^s + o-)e—(2p3-2pi) ^ 

= -{K\ + 0")e-"(2p^-2Po) ^ 


where O'' := 0{e-''^). 
Setting 


= Po + Po + Po > 0 , (Tp, = Pi + Pi + Pi , 

the leading order behaviour of the metric components, and those of its inverse, is 


Poo 

= - ( 

g-2<7pc 

(1 + 0 ^) ^ 

0 , 

POi 

= 0, 




Pll 

= e~ 

-2piT- 

2 P^(l+0"), 


P22 

= e~ 

-2p^r- 



P33 

= e~ 

-2Po'r- 

2^°(1 + O")^0, 


P12 

= e~ 

-2p^g 


0 , 

Pl 3 

= e~ 



0 , 

P23 

= e~ 

-2Po-e 


0 , 


^00 ^ _^2a^.r+2ap,^-^j^QU^^ 

5°* = 0, 

gii = e"P°"+2^°(l + O")^0, 

= e^P>+2Pl ^1^0''), 

f3 = e2Po^+2/3c(l + 0'^), 

gl2 ^ _g2pirg2/3^ (^1^ + ^ 0 , 

gi3 = +0") ^ 0, 

f3 = _e2plr^2p^j^2^Q.^^ 


(E.8) 


(E.9) 


where all decaying ones have been marked with “—>■ 0”. 

The Christoffel symbols are defined as 

fr /37 “ 2 ^ i9o-'y,P + gP(T,j ~ P/ 37 ,O') , (E.IO) 

and show the following behaviour: 

r %0 = -'^Po(l + 0{e-n), rO,o = + a^„,,)(l + 0 (e—)), 

TV = -Po(l + 0(e-"")) (no sum), = -pyOo+0o)e^(.pl+plP(l + 0(e-"")) , 

r °22 = -p2e2(/3o+/3o")e2(p)+Po)-(i + 0 (e-"")), r 033 = -p3e2(/3)+/3o")e2(p)+Po)-(i + 0 (e-"")), 
r°2i = (C + 0'^)e2fr°+P°)^ , r°3i = {C + O^)e2fro+Po)^ ^ 
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r032 = (C + 0'^)e2(Po+Po)^ , 

r\i = (C' + o>, 

T\^ = {C + OnT, 

rigo = (C + 0'^)e2(Po-p?)^ ^ 0, 
r^32 = —)■ 0 , 

= (C + 0'^)re-2(P°+P°)^ ^ 0 
= (C + 0'^)re2(-P°+P°)^ 
r 222 = (c + o")r, 
r233 = (C + 0'^)r, 
r233 = (C + 0")Te2(P°-Po)") ^ 0, 

r\o = (c + o'^), 

T\^ = {C + 0^), 

= (C + 0 '^)re 2 (-P°+P°)^ , 

T\^ = {C + Onr, 


r^oo = (C + 0^)re-2(Po+Po)^ ^ o , 
T\o = ^ 0, 

T\2 = iC + 0^)Te2(P°-P°)^ ^ 0 , 

T\, = {C + Onr, 

ri33 = {C + 0^)re2(P°-Po)^ ^ 0 , 

T\, = {C + On, 

T\ = {C + Onr, 

r^32 = (c + onr, 

T\o = (C + 0")Te2(-P°+P°)" ^ 0 , 
= (c + o^)Te2(-p°+p°)^ 

= {C + 0^)Te2(-P°+P°)^ , 

T\, = iC + Onr, 

T\, = {C + Onr, 


with C denoting r-independent quantities, possibly different for different components. For most 
components an explicit first order term is given above. This is the case if the highest order term 
appearing in (E.lOl does not vanish and therefore only the first order terms of the metric, given 
in (E.9) contribute. For the other terms, e.g. r^ 20 : the highest order term cancels and lower order 
terms of the metric could potentially become important, necessitating a more detailed analysis. 
Here the behaviour of these terms is simply given as decaying faster than the vanishing highest 
order term, as this is sufficient to determine the behaviour of the Kretschmann scalar in the 
following. 

From the Christoffels we calculate the behaviour of the coordinate components of the Riemann 
tensor, defined as 


— ^acr (F /3(5,7 


r tT I T^cr T^e T^cr -pe \ 

136 — ^ 6e^ /37j • 
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Their behaviour 

is given by 









.Rioio = 


0’')e- 

t2Po 

5 


.R 2010 = 

--{C 


0’')e- 

r 2 pl 

—)■ 0 , 

-^2020 = 

(C^ + 

0’')e- 

^ 0 , 

-R 2110 = 

--{C 


0'')Te 

— r2p 

1 

0 

5 

.^2120 = 

(C-T 

C»^)re 

-r 2 pl 

—)■ 0, 

.R 2121 = 

-{C 




-^3010 = 


0")e- 

^ 0 , 

.R 3020 = 

--{C 


0")e- 

r 2 pl 

—)■ 0, 

-^3021 = 

o(Te~ 


0, 


-^3030 = 

--{C 


0’')e- 

r 2 pl 

—)■ 0, 

-^3110 = 


C»^)re 

-'r 2 pl 


.R 312 O = 

-{C 


0^)re 

— r 2 p 

2 

° 0, 

-^3121 = 

{C + 



-^3130 = 

= {C 


0^)Te 

— r 2 p 

3 

° 0, 

-^3131 = 




.R 32 IO = 

--{C 


0’')Te 

— t2p 

2 

° 0, 

-^3220 = 

{C + 

0^)Te 


—^ 0, 

-R 322 I = 

--{C 



-^3230 = 


0 ^)re 

-r 2 pl 

—)■ 0 , 

.R 323 I = 

-{C 


^0, 

R3232 = 


0 ’'yr 2 pl 

0 . 








The Kretschmann scalar, defined as RajS-ysR’^^^^, behaves as 

^t4{pI+pI+pI) 


with 


K = = {C^ + 0(e 

16g4o-/3o 


l)e- 


= 


{pipin{pif+pipi+{pi?^ 


{pI+pIY 


(E.ll) 

(E.12) 


As the inequalities (8.14) require p^> pi > Q the coefficient is positive and the Kretschmann 
scalar diverges as r —>■ oo for all constructed solutions. 


Appendix F. Solutions including a cosmological constant 

Solutions of the form presented above can also be constructed with a nonzero cosmological 
constant. Here we will show the resulting changes in the arguments above. 

The presence of a cosmological constant changes the action to 


S[g^Y = / d^xy^(i? - 2A). 


(F.l) 


This leads to an additional term 2Kg in the Hamiltonian. As the determinant g of the spatial 
metric behaves asymptotically as 

(F.2) 

the term decays exponentially. It contributes an additional term to the evolution equation for 


TTa, (8.2b). This equation includes the diverging prefactor but the conditions (8.9) guarantee 


that the new term decays fast enough to compensate it. Therefore the presence of a cosmological 
constant introduces no new conditions on the free functions from the evolution equations. 

The cosmological constant also appears in the Hamiltonian constraint, but not in the mo¬ 
mentum constraint. As the additional term in the Hamiltonian decays it does not change the 
asymptotic Hamiltonian constraint or the condition on the pi arising from it. 
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Finally, A appears in the derivation of the evolution equations for the constraints in Append!^ 
[d| Equations (D.8) and (D.IO) change respectively to 


i? = -^Goo + 2gA + 0(7V'=), 
+ 0{N'^) . 


(F.3) 

(F.4) 


The additional terms containing A cancel in the final evolution equations (D.ll) and (D.12). 


E 


as 


Therefore the arguments regarding the constraints in section 7.3 remain unchanged. 

We conclude that the solutions with the asymptotic behaviour given in [Appendix 
obtained from the Fuchs theorem above, exist for all values of A, with the free functions and the 
associated conditions unchanged from the A = 0 case. 
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